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Abstract. A continuous projective measurement of a quantum system often leads 
to a suppression of the dynamics, known as the Zeno effect. Alternatively, generalized 
nonprojective, so-called "weak" measurements can be carried out. Such a measurement 
is parameterized by its strength parameter that can interpolate continuously between 
the ideal strong measurement with no dynamicsthe strict Zeno effect, and a weak 
measurement characterized by almost free dynamics but blurry observations. Here 
we analyze the stochastic properties of this uncertainty component in the resulting 
observation trajectory. The observation uncertainty results from intrinsic quantum 
uncertainty, the effect of measurement on the system (backaction) and detector 
noise. It is convenient to separate the latter, system-independent contribution from 
the system-dependent uncertainty, and this paper shows how to accomplish this 
separation. The system-dependent uncertainty is found in terms of a quasi-probability, 
which, despite its weaker properties, is shown to satisfy a weak positivity condition. 
We discuss the basic properties of this quasi-probability with special emphasis on 
its time correlation functions as well as their relationship to the full correlation 
functions along the observation trajectory, and illustrate our general results with simple 
examples. We demonstrate a violation of classical macrorealism using the fourth-order 
time correlation functions with respect to the quasi-probability in the twolevel system. 
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1. Introduction 

The continuous projective von Neumann quantum measurement |TJ leads to a 
suppression of the dynamics, known as the quantum Zeno effect (QZE) [2]. To 
escape this problem, modern quantum measurement theory offers generalizations of the 
projective measurement to so-called positive operator- valued measures (POVM) [3111], 
where a projection is replaced by a softer operation such as a Kraus operator [5]. Such 
operators can describe not only projective measurements but also weak measurement, 
in which case the action of the POVM leaves the state almost unchanged. By virtue 
of the Naimark theorem [B] POVMs are equivalent to projective measurements in an 
extended Hilbert space that includes additional detector degrees of freedom. The effect 
of a continuous application of Kraus operators, which correspond to a time-continuous 
measurement can be described by stochastic evolution equations such as Lindblad-type 
equations [7] for the system density matrix or Langevin equations for individual system 
trajectories, physically describing irreversible effects like decoherence and decay affected 
by the measurement process. 

Weak measurements [H] make it possible to escape the QZE by paying a price in 
terms of an imperfect detection. In the extreme case the dynamics of the system is 
(almost) free but the measurement outcome is obscured by a large detection noise. 
This is similar to the problem of a quantum linear amplifier, which can amplify both 
complementary noncommuting observables, like x and p, but only if accompanied by 
a large noise [9J. The interpretation of weak measurements of correlation functions is 
sometimes paradoxical: one must either accept unusually large values of the physical 
quantity [8 J or replace the probability by a quasiprobability [TP] . Weak measurements 
are also very useful in quantum feedback protocols [IT] . 

The QZE lies at the strong limit of a spectrum of measurements whose strong/weak 
character can be changed continuously p21 [131 EH HSJ HE] UT] HE] EE], using e.g. a 
Gaussian POVM [201 EI]- The Gaussian POVM is also the key element of the continuous 
collapse interpretation of quantum mechanics [22] • These models lead to various types 
of expressions for time correlation functions [23]. Here, for the first time we use 
Gaussian POVM for continuous measurement to describe higher order symmetrized time 
correlation functions. Such functions are known for the two-time case of a harmonic 
system |20j or in the weak measurement limit for the two-level system [21] (but not 
in general) and are necessary to explain the recent experiment that shows nonclassical 
behavior of time correlation functions in a two-level system [25J. The calculations are 
faciliated by making a deconvolution of the outcome time trace probability into the 
probability component associated with the white detection noise and a quasiprobability 
that describes the intrinsic system uncertainty. Such a deconvolution has the advantage 
that we can make use of basic properties of the quasiprobability, e.g., the weak positivity 
[26J, which states that the second order correlation function matrix is positive definite. 
Our scheme provides a unified and concise treatment of weak measurements and the 
QZE, pointing out the general trade-off between measurement and decoherence. By 
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comparing the average signal to the associated noise we also establish limits on the 
uncertainty of the outcome and its dependence on measurement characteristics. 

The time correlation functions obtained by our approach provide a convenient route 
for the analysis of uncertainty properties of systems undergoing weak measurements. 
Taking a two- level system as an example, a single nondemolishing measurement of 
an observable not commuting with the Hamiltonian is not possible in both time and 
frequency domain although the latter gives a better signal-to-noise ratio. Although this 
is intuitively clear, using our approach, it is possible to establish and compare bounds 
on the outcome uncertainty. For another simple example, the continuous position 
measurement of a harmonic oscillator, we show that the same measurement procedure 
does not lead to the QZE. Instead, the continuous measurement leads to unbounded 
growth in noise, in agreement with the general observation [2T] and in analogy to anti- 
Zeno effect [27] . 

The proposed separation has another important consequence. If we assume classical 
macrorealism in quantum mechanics then the statistics of the outcomes with the 
detection noise subtracted in the limit of noninvasive measurement should correspond 
to a positive definite probability. In contrast, we show that macrorealism assumption 
is violated by demonstrating that our quasiprobability is somewhere negative. Such 
violation has been recently demonstrated experimentally j2Sj. In fact, if we additionally 
assume dichotomy or boundedness of the quantum outcomes, the violation can occur 
already on the level of second order correlations of a single observable as shown by 
Leggett, Garg and others [2UEH], but also indirectly, subtracting the unavoidable (and 
necessarily divergent) noise [29]. However, as follows from weak positivity, without 
these additional assumptions, second order correlations are not sufficient to violate 
macrorealism. Instead, one needs at least fourth order averages to see this violation. In 
this paper we demonstrate, that a special fourth order correlation function in the two- 
level system, reminiscent of the Leggett-Garg proposal[28j. can reveal the negativity of 
the quasiprobability in this case and consequently can be used to violate macrorealism, 
without having to make any additional assumptions. 

The paper is organized as follows. We first define the continuous Gaussian POVM 
and obtain the probability distribution for the continuously measured observable. 
We then make the deconvolution of this probability into detection noise and a 
quasiprobability and introduce a formalism for evaluating time correlation functions. 
With that we are able to prove the weak positivity. We show that the time evolution 
associated with the quasiprobability can be formulated either as a quantum Langevin 
equation driven by a white Gaussian noise, or a Lindblad-type master equation for the 
non-selective system density matrix. We also show how the required time correlation 
functions can be calculated from these stochastic equations. Next, we demonstrate 
the general trade-off between dynamics and measurement, taking a two- level system as 
example, and discuss the behavior of the average signal and the noise in these prototype 
system. Then we construct the Leggett-Garg inequality without assuming dichotomy 
or boundedness of the measurement variable. Finally, for completeness, we discuss the 
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harmonic oscillator case and show how and when the Zeno effect emerges within our 
formalism. Several instructive proofs of formulas are presented in Appendices. 

2. Quasiprobability and weak positivity 

We begin by introducing a general scheme of continuous measurement and describe its 
properties. For a given system characterized by a Hamiltonian H and an initial system 
state given by a density matrix p, we consider the measurement of one, generally time- 
dependent, observable A. A description amenable to continuous interpolation between 
hard and soft measurement can be formulated in terms of Kraus operators [5j |8]. We 
assume a Gaussian form of the Kraus operators, whereupon the state of the system 
following a single instantaneous measurement is given by 

Pl (a) = K(a)pk(a), (1) 

K(a) = (2\/n) 1 / A e-^ a - A ^. (2) 

Note that in Q the non-negative definite operators p and p\ represent the states 
of the system just before and just after the measurement. The probability that the 
measurement of A gives the outcome a is given by [3] 

P(a) = Trpi(a), (3) 

which is normalized, J da P(a) = 1. The Kraus operator (|2j depends on the parameter 
A, which characterizes the weakness of the measurement. For A — > oo, we recover 
a strong, projective measurement with an exact result but a complete destruction of 
coherence, while A — > corresponds to a weak measurement with almost no influence on 
the state of the system, pi(a) ~ p, but very large measurement uncertainty of the order 
of ~ 1/A. The probability distribution ([3j is consistent with projective measurement 
scheme, namely (a) = TrAp. 

Let us imagine a continuous sequence of meters interacts with the system. The 
meters are prepared with a Gaussian wave function, the interaction is proportional to 
the product of the system observable A and the meter momentum, and the position of 
each meter is read out after the interaction The post-interaction position of the meters 
is the measurement result et(t)[U HH |23l 120] . 

Repeated measurements of this type can be described by applying such Kraus 
operators sequentially, separated by time steps At. In what follows we make the 
reasonable assumption that for a given measuring device ("meter") the weakness 
parameter A is inversely proportional to the measurements frequency, i.e., 

A = XAt (4) 

with constant A. In the continuum limit, A, At — > 0, we obtain (Appendix A) the Kraus 
operator as a functional of a{t) 

K h [a(t)} = e^ At K [a (t)] = CTe' * MQ-Mt))** (5) 
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where a(t) is the measurement outcome, A(t) is the operator A in the Heisenberg rep- 
resentation with respect to the Hamiltonian H, A(t) = exp (iHt/hj Aexp (^—iHt/hj, 
Tdenotes time ordering (later times on the left) and C is a normalization factor. Note 
that Kh[o] is the Heisenberg representation of K[a). The analog of (J3J) is the functional 
probability 

P[a]=Trk*[a]k[a]p (6) 

which satisfies the normalization J DaP[a] = 1. Whenever some functional measure D 
is introduced here, we tacitly include all proper normalization factors in it. 
It is convenient to write (j5l) as a Fourier transform 

K h [a\ = J ^ Te I^H{t)(Mt)-a{t))-^{t)/AX] (7) 

so that 

P[a] = Tr [k\a]pk\a]) = Tr (J D +re / d *^+W( i «- a W)-^«/ 4A ] x 

p J D0„Te^ d *^- (i)(i( * ) " a{ * )) " 9i - ( ' )/4A1 ^ (8) 

where Tdenotes inverse time ordering (later times on the right). Changing integration 
variables according to x — 4 ) + + 0- an d 4> = (<t>+ — 0-)/2 we can write 



P[a] = J D 0e-/ d ^«/ 2A J B X e-^ dt ^ 8X e-^^ a ^ dt (9) 
x Tr j- e IKx(t)/2+Ht))Mt)dt p f e Ji( x (t)/2-Ht))A(t)dt 

The last line can be written alternatively as (see Appendix B) 

Tr Te^ xit)A(t,{t)dt/2 p fe i f x{t)A ' t ' (t)dt/2 (10) 

where A$(t) denotes the operator A in a modified Heisenberg picture, namely with 
respect to the Hamiltonian H — h<p(t)A. 

([9]) and (10) describe the outcome of the continuous measuring process in terms of 



the probability distribution functional P[a(t)] of the observation function a(t). This 
distribution reflects the quantum uncertainty, the modified system time evolution caused 
by the measurement (the backaction effect) and the uncertainty associated with the weak 
measurement that can be thought of as reflecting detector noise. A more transparent 
view of these contributions is obtained by separating the latter, system independent 
contribution from the system dependent effects. This is achieved by considering the 
moment generating functional (MGF) M[x] = e s ^ x \ where S [%] is the cumulant 
generating functional (CGF), given by 



M[ X ] = e 5M = J Da x(t)a(t) P[a) = (11) 

e -/dt x 2 (t)/8A f D e -/d^ 2 (t)/2A Tr j-jfx(t)Ai(t)dt/2 p j- e i/x(t)^(*)dt/2_ 
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The CGF can be divided into two parts S[x] = Sd[x] + Sq[x] with 

S d [ X ] = - j dt X 2 (t)/8X (12) 

and 

e S q [x] _ f^ e -fdt4>Ht)/2X Tl j- e if x (t)A^t)dt/2pf e iJx(t)A4t)dt/2^ /-^ 



Note that S[0] = S d [0] = S q [0] = 0. On the level of probabilities this decomposition 
corresponds to the convolution 



P[a] = J Da'P d [a-a'}P q [a'], (14) 

where 

P d [ a ] = [ B X ef dtMt)a{t)/i - xHt)/sx) oc e - 2A ^ 2( ' )d * (15) 



corresponds to a Gaussian noise with zero average and correlation (a(t)a(t')) d = 
5(t — t')/4X that may be interpreted as the noise associated with the detector, and 
where 

P q [a] = J B X e-^ xit)ait)dt e s " [x] (16) 

is a distribution associated with the intrinsic system uncertainty as well as the 
measurement backaction. It is normalized, J DaP q [a] = 1, but not necessarily positive, 
and will be referred to as a quasi-probability [TO], [301 El] - m the limit of weak, noninvasive 
measurement, A — > 0, P d diverges while P q has a well defined limit 

e S g [x] ^ Tr j- e ifx(t)A(t)dt/2 p j- e if X (t)A(t)dt/2_ 

Consider now this distribution P q . First note that while it is not a real probability 
functional, it is possible to calculate moments () and cumulants (()) with respect to this 
measure as partial derivatives of the quasi-CGF, respectively 

S n expS g [ X } 



{afa) ■ ■ -a(t n )} q = 
((a(ti) ■ ■ -a(t n ))) q 



Six (ti) ■■■Six (tr 

s n s q [ x ] 



(18) 



Six (ti)---tfix (t n ) 
In particular, for t n > . . . > t 2 > t±, 

(a(t)) q = [ D A 0Tr[i^(t)p], (19a) 



(a(t 1 )a(t 2 )) q = J D A 0Tr[{A (ti),^(t 2 )}p]/2, (196) 
(a(t 1 )a(t 2 )a(t 3 )) q = J D^Tr^i^), {^(t 2 ), i (t 3 )}}p]/4, (19c) 
(a(h) ■ ■ ■ a{t n )) q = [ D^TriiA^h), {{A^t 2 ), ■ ■ ■ A^(t n )} ■ ■ -}p]/2"-\ 
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(see Appendix C), where we have defined = D0e J d * * 2 ( t )/ 2A . Here and below we 
use standard notation {A, B} = AB + BA and [A, B) = AB - BA. 



Secondly, from Eq. (196) follows the important so-called weak positivity property 



of second order correlations |26] 

(20) 



(F 2 [a]) q = J D0e-^ 2 W/ 2A TrF 2 [A ]p > 



for F[a] = J dt(f(t)a(t) + g(t)) and arbitrary functions / and g. It can be interpreted as 
a generalization of the Robertson-Schrodinger uncertainty principle [32J . This property 
has an important implication that no test based solely on maximally second order 
correlations can reveal the negativity of the quasiprobability. First and second order 
correlations can be represented by a completely classical, positive Gaussian probability 
distribution 

P' q [a] oc exp f- j dtdt'5a(t)f- l (t,f)Sa(t r )/2 S j , (21) 

where 5a(t) = a{t) — (a(t)) q , f(t,t') = (5a(t)5a(t')) q and / _1 is its inverse defined 
by f dtf(t', t)/ _1 (t, t") = 5(t' — t"). The weak positivity guarantees that both / and 
f^ 1 are positive definite and consequently P' q is a correct real probability distribution. 
To check that P q differs from P' and demonstrate its negativity one needs higher order 
correlations or additional assumptions (e.g. boundedness or dichotomy of a as it happens 
in Leggett-Garg inequality [28J). 

To end this Section, we consider the special case in which the Hamiltonian 
commutes with A (or the noncommuting part is negligible during the interesting 
timescale). Furthermore, let us take the initial state of the system to be an eigenstate 
| a) of A, i.e. p (t — 0) = |a) (a| , A\a) = a\a). Consider a measurement performed during 
the time interval to, 

a = (l/t ) f ° dt a(t) (22) 



In this case we find (Appendix D) that (a) = (a) q = a and ((5d) 2 ) = {{5a) 2 ) d = l/4A£o 
with 5X = X — (X). We can see the intuitively expected effect of an increasing 
measurement duration to lead to an improved signal to noise ratio with time, which 
goes as 

(fi> 



2a v / At - (23) 



Thus, even the weakest measurement (small A) turns into strong a one, if performed 
often enough for sufficiently long time. 

3. Representation by stochastic evolution equations 

Turning back to the general case, we note first that the correlation functions associated 



with the quasiprobability P q [a(t)], given by Eqs. (19a 19c), can be calculated from the 
Heisenberg equations 

dB+(t)/dt = (i/h)[H+(t) - <j)(t)Mt), B+(t)] > ( 24 ) 
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where A represents the measured variable while B is any system operator. In particular 
dA^/dt = (i/h)[H^t),A^t)l 

dH^/dt = </>(t)dA4,(t)/dt. (25) 
We can solve these equations for a general stochastic trajectory <f>(t), then take the 



averages as defined by Eqs. (19a 19c), over a Gaussian distribution of such trajectories. 
The correlation functions obtained in this way coincide with the ones derived directly 
from the CGF Eq. (16). If H = p 2 /2m + V(x), with [x,p\ = ihl, and A = x is the 



position operator, the Heisenberg equations for x<j> and are 
H^=pl/2m + V{x 4> ), 

^ = fc/m, (26) 



(26) is a quantum Langevin equation in which the quantum dynamics is augmented by 
a zero centered white Gaussian noise, (4>(t)) = 0, (4>(t)(j)(t')) = \5(t — £'). Closed form 
solutions of this equation can be obtained for the harmonic oscillator, a case we discuss 
below. 

Alternatively, the stochastic dynamics affected by the continuous measurement 
process may be described by a Lindblad-type master equation [7j for the non-selective 
system density matrix. The latter is defined by 

/■a(t) pa(t) 

p(t) = / Dap[a] = / DaK[a]pK j [a] (27) 

Ja(0) Ja(0) 

where the integral is over all observation trajectories between times and t. It is shown 
(appendix E) to evolve according to (using a Liouville-superoperator L) 

^ = Lp := [Hj]/ih-X[A, [Aj}}/2. (28) 

In the representation of eigenstates of A, 

P = ^2paa>\a)(a'\, (29) 

a,a' 

Lp a ,a> = 7T ^ {H ab p ba , - p ab H ba <) - A (a - a') 2 p aa >, (30) 

b 

showing, as is well known [lj and as may be intuitively expected, that the measurement 
damps the off-diagonal terms (a ^ a') with the rate proportional to the measurement 
strength. Note that if some eigenvalues a are degenerate, then the corresponding off- 
diagonal elements of p are not damped. 

Together with the Liouville-Lindblad superoperator L we define the corresponding 
evolution superoperator U(a,b) = Texp J? Ldt. It can be then shown (Appendix F) 
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that the correlation functions (19a 19c) are given by 

(a(t)) q 



Tr 

(a(ti)a(t 2 ))q = Tr 
(a(ti)a(i 2 )a(*3)) g 



AU(t,0)p 

Au(t 2 .t t ) Ant 

■■ Tr 



(a(ti) • ■■a(t n )) c 



Tr 



AU(hM)AU(t 2 M)AU(nQ)~p 
AlJ(t n , t n -x) ■ ■ ■ AU(t 2 , t x )AU{t x , 0)p 



(31a) 

(316) 
(31c) 



where AB = {A, B}/2. Note that in (|31af31c|), p = p (t = 0) = p. (j31a§31c|) provide a 
more convenient route for the evaluation of these correlation functions. 

In the following Sections we apply this general formalism to the two simplest 
quantum systems, the two-level system and the harmonic oscillator. 

4. The two-level system 

Consider a two-level system defined by the Hamiltonian 

H = Hu& x /2 (32) 
and suppose that the system is in the initial state 

p(t = 0) = (i + <7 z )/2, (33) 

where a denotes Pauli matrices and 1 is the corresponding unit operator. Left 
uninterrupted, the system will oscillate between the two eigenstates of a z , a process 
analogous to Rabi oscillations in a harmonically driven system. We focus on the 
measurement of A = d z and denote the measurement outcome by a(t) = a z it). We 
pose the following questions: Can the oscillatory time trace of a z be observed? How 
does the measurement process affect this oscillation? Is the oscillation visible in a single 
run of an experiment or only as a statistical effect - average over many runs or many 
copies of the same experiment? The latter question is particularly relevant in light of 
the growth of activity in single molecule spectroscopy. 



To answer these questions we start by writing the action of L, (28), in the basis of 
Hermitian operators (a x , & y , a z ). In a compact notation it reads, 

L(x& x + ycr y + za z ) = u(ya z - za y ) - 2\(x& x + ya y ) 

and LI = 0. Next, expressing the operation of U (t, 0) on b z by 

cr z (t) = U (t, 0)a z = x{t)a x + y(t)a y + z(t)a z 

-2Xx; dy/dt = 



and using (34) and (J35J) in (28) we find dx/dt - 
dz/dt = coy, which, for z(t = 0) = 1, x(0) = y(0) = yields 



(34) 
(35) 

(cuz + 2Xy), and 



z(t) = e~ xt {cos(ttt) + Asin(nt)/n], y(t) = dz/dt, x(t) = 0, 



(36) 



where = y/tu 2 — A 2 . This allows to write down the relevant averages (see Appendix 
G), namely 



Mt)) q = z{t), (a z {t)a z {t')) q = z{\t - t'\). 



(37) 
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The last line is known in the existing literature only for the stationary case (t, t' — > oo) 
in the weak measurement limit A — > |23J[T7]. It is interesting to note that although the 
system under consideration is not in a stationary state and in fact evolves irreversibly, 
this correlation function depends only on the time difference t' — t and remains finite 
when this difference is constant while both t and t' increase. 

Recall that ( [573] ) and flDl| imply that (a z (t)) = (a z (t)) q , while (a z (t)a z (tf)) = 
(a z (t)a z (t')) d + (a z (t)a z (t')) g , which implies (5a z {t)5a z {t')) = (5a z {t)5a z {t')) d + 
{Sa z (t)Sa z (t')) q . In the limit ^ < w we see (cf. (36) and (37)) clear oscillation of 

(o~ z {t)). However, in a single run this signal cannot be distinguished from the noise. 

-l 



Indeed, defining as in (22) o z 
Therefore 



(1/to) L° dta z (t), we obviously need to take t <w 



((5a z ) 2 ) > (((5a z ) 2 ) d = 1/toX » 1 (38) 

The large detection noise covers the signal. This implies that Rabi oscillations cannot 
be seen in a single run. 

The above result was obtained in the time domain. We can also ask whether the 
Rabi oscillation is visible in the frequency domain. This would imply seeing the peak 
in the Fourier transform 



Or 



to 



0) = (2/t ) / dtcos(pt)a z (t), 



(39) 



where to is a time much larger than the oscillation period, but obviously much smaller 
than the damping time: u^ 1 <C to <^ A -1 . From (36) and (37), the peak intensity 
is (a z [y = Q)) ~ 1. On the other hand, under the measurement conditions the white 
detector noise satisfies 

((5a z (n)) 2 ) > ((5a z (n)) 2 ) d = l/2t Q A » 1 (40) 

implying that, again, the noise exceeds the signal and a peak in the frequency domain 
will not be seen. This time, however, the signal-noise ratio is not as bad as in the time 
domain because t can be longer. 

We conclude that Rabi oscillation cannot be seen in a single run/copy of the 
experiment but only in a statistical average. The sample size, that is the number 
of runs/copies needed for this average is of the order (t A) _1 , where A -1 ^> t ^> u^ 1 
in the frequency domain and to ou~ l A -1 in the time domain. In the overdamped 
regime, A ^> u, one can see the quantum Zeno effect, discussed below in Section 7. 



5. Leggett-Garg-type inequality 

The limit A — > is consistent with the noninvasive measurement because the backaction 
vanishes. In this case the negative quasiprobability demonstrates the violation of 
macrorealism even for a single observable, as shown by Leggett and Garg [2H]- In 
violations of this type, it is essential to subtract the large detection noise, whose 
uncertainty must always diverge and prevent any real violation [29]. The common 
confusion about the noninvasiveness condition is caused by the fact that two-time 
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correlations are numerically identical for the quasiprobability in the limit A — >■ and 
the instant projections (invasive because of collapse) for initial p ~ 1. The equality 
still holds in the case of many times if the observable satisfies A 2 oc 1. The analysis 



above has used second order correlations that, as stated in (20), are not sensitive to the 
quasiprobabilistic nature of the distribution. The violation of the well known Leggett- 
Garg inequality [28] needs only second order correlations but requires the additional 
assumption of bounded observables which is effectively equivalent to higher order 
correlations (e.g. the dichotomy A = ±1 is equivalent to measuring ((A 2 — l) 2 ) = 
which requires the fourth-order correlator (A A )). Without this assumption, the 
quasiprobabilistic nature is however revealed in higher order correlations. To see this 
we take p = 1/2 and consider the following quantity 

X[a] = a z {0)a z {7r/u) + a z {-7r /2uo)a z {-K /2uo) + 2. (41) 

The fourth order correlation (X 2 ) q is given by 

(X 2 ) q = 6 + e~ A7r/a; [l/r 2 + (10 - 1/r 2 ) cos(Trr) + 10A sin(7rr)/wr], (42) 

where r = a/1 — (\/u) 2 . The behavior of (X 2 ) q is shown in figure [TJ In the limit 
of strong measurement (X 2 ) q = 16. The origin is the QZE - the evolution is frozen 
by the measurement and so c z (t) does not depend on time, which results in X = 4. 
In the opposite limit of noninvasive measurement (X 2 ) q = —2 and it crosses zero at 
u/X ~ 11. This implies that for sufficiently small A the classical inequality (X 2 ) q > 
is violated so the function P q is not positive definite and as such cannot describe a 
usual probability. Note, however, that (a) it contains the relevant physical information, 
discarding the irrelevant detection noise; (b) by itself, it cannot be directly measured, 
namely correlations such as (X 2 ) q are not directly measurable since the real probability is 



the convolution (14); and (c) the actual detected observable certainly satisfies (X 2 ) > 0. 
However, an independent determination of the detection noise should be experimentally 
feasible and allows to find the negativity of (X 2 ) q after the noise has been subtracted. 

6. The harmonic oscillator 

For completeness, we also consider another much studied simple problem - continuous 
position measurement, A = x, in a system comprising one harmonic oscillator 



described by the Hamiltonian H = p 2 /2m + mu 2 x 2 /2. Eqs. (26) become 

dx^/dt = p^/m, 

dp^/dt = -mu 2 ^ + h(f)(t), 

where (f) represents the zero-centered white Gaussian noise, {4>(t)(j)(t')) = \5(t — t'). We 
note in passing that this quantum Langevin equation yields the Fokker-Planck equation 
for the Wigner function [33] 

W(x,p)= [ p^ e -^ x -' 1XP Tipe^ x+ixP (44) 
J (2vr) 2 
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CM 

X 




Figure 1. Demonstration of the violation of the Leggett-Garg-type inequality 
(X 2 ) q > as a function of measurement strength A. The function starts from the 
value 16 for A 3> uj as expected for the QZE, crosses the classical bound and tends 
to —2 in the limit of weak measurement A <C w. 



in the form 



dW(x,p,t) 2 dW 
- = mu x 



p dW Xh 2 d 2 W 
+ 



(45) 



dt dp m dx 2 dp 2 

However, in what follows we calculate directly the required correlation functions. Solving 



(43) we get 



xAt) = x(0) cosut H sinut+ / smwtt — t ) n<p(t), 

mu; J muj 

ft 



(46) 



Ptj,(t) = p(0) cosut — mojx(0) sinwt + / dt' cosu(t — t') hcj)(t'), 

Jo 

This implies that the oscillations of the average position are undamped, 

(x (t)) = (x (0)) cos (ut) + (my) -1 (p (0)) sin (tut) 

independently of the detection strength. 

Turning to the noise term we first note that, as before, the detector noise combines 



(47) 



additively with the the correlation functions obtained from (46). The latter take the 
form 

(Sx (t) 5x (t')) q = (Sx (t) 5x (O) + fx (t, t') (48) 

where (8x(t)5x(t'))a is the free correlation function obtained in the limit A — > or 
equivalently </> — > 0, that is, by ignoring the last (noise) terms on the RHS of ( |46| ), 

(6x(t)6x(t')) = (5x(0)Sx(0)) coswtcoswt' + (49) 
(5x(0)5p(0)) H /(mu;)" 1 sinw(t + t') + (5p(0)5p(0))(muj)~ 2 sin sin ut' 
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with the Wigner-ordered average (2xp)w 



correlation function associated with the noise terms in (46) 
/a(M') = 



Tr[p{£,p}], and where f\(t,t') is the 

(50) 



Xh 2 



-[min(t,t') cosw(t - t') + (sina;|t - t'\ - sinu;(t + t'))/2u}. 



2(mu) 2 

This measurement-induced correlation function represents the backaction effect of the 
measuring process. It depends on the detector strength and the parameters of the 
dynamics but not on the initial state of the oscillator. Moreover, because of the Gaussian 
nature of (f>, it contributes solely to the second cumulant ({x(t)x(t'))) , leaving all the 
other unaffected. As expected, it vanishes in the limit A — > 0. However, the most 
striking feature in (50) is the growth of noise with time, as expressed by the first term 
in @. 

In analogy to the two-level system we discuss the behavior of the short time 
(t <C average x and the long time (t ^> uj^ 1 ) Fourier transform x, defined by 

the analogs of (22) and (39), respectively. In both limits we are now free to choose to\ 
because, in contrast to the 2-level case, the averaged oscillations, (47), are not damped. 

Consider first the time-domain observation. For t A <C 1 the uncertainty of x 
is determined by the detection noise as in (38), while in the opposite limit it will be 
dominated by the backaction (50). In either case the noise exceeds the signal. 

In the frequency domain, for x(u) = 2 J *° dtcos(ojt)x(t)/t , we get using (47 50), 
the peak signal 

(x(w)) = <x(0)), (51) 

and the intrinsic and backaction noise components 

Xt h 2 



((Sx(co)) 2 ) a ^((5x(0)) 2 ) + 



(52) 



The total uncertainty originated 



(53) 



6 (muS) 

to which we need to add the detector noise (2t A) -1 . 
from the detector satisfies 

Xt h 2 1 k 

u | > 

6 (mw) 2t A y/Smui 
with the lower bound (obtained as minimum of the LHS with respect to t X) independent 
of A and t . Obviously (x (0)) can be chosen large enough for the signal to dominate the 
noise at intermediate times, but the noise will always exceed the signal at long enough 
times. As always, the signal-to-noise ratio can be improved by repeated measurements. 

The fact that the backaction contribution ( 50 ) to the noise grows with time reflects 
the continuous pumping of energy to the system affected by the measurement process 
[55] . This does not happen in the two-level system because of its bounded spectrum, 
still also in that system the temperature grow to infinity {pit) — >■ (1/2)1) as implied 
by Eqs. (35), (36) and (37). This unlimited growth can be avoided by assuming that 
the measurement process also involves some friction [211 E!^ 



Indeed, measurement, 
even classical, means extraction of information out of the system, so that without 
compensating friction its entropy must increase and so does the temperature. 
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7. The quantum Zeno effect (QZE) 

For completeness, we show now how the QZE emerges within the present formalism. So 
far we have focused on weak measurements, represented by small A. The opposite limit, 
A — > oo, represents the strong measurement case. In systems characterized by a single 
timescale w -1 , strong and weak measurements are quantified by the inequalities A 3> cu 
and A <w, respectively. 

Consider the two-level system discussed in Section 4. For A > u its dynamics is 



given by the overdamped analog of Eq. (36), Vt = i\/\ 2 — to 2 . In the extreme strong 



measurement case, A ^> uj, z{t) ~ e~ u t / 2X and the decay slows down as A — > oo[Hl 
This corresponds to the QZE where the system is almost frozen by the measurement 
reaching its equilibrium state z = only on the timescale t ~ X/u> 2 . 



For a position measurement in the harmonic oscillator case, we have seen, Eq. (47), 
that the average position oscillates regardless the strength of the measurement. This 
implies that the Zeno affect is absent in this system, as is well known [T4] . On the other 



hand, for any measurement strength, the detector- induced backaction noise, Eq. (50), 
increases without bound at long times at a rate that increases with A. Already for short 
times we get f\(t,t) ~ Xh 2 /3m 2 , and backaction adds fast diffusion in the phase space. 
This is somewhat analogous to the anti-Zeno effect |27j . 



8. Conclusions 



Gaussian POVMs, here represented by Kraus operators, were used in this paper to 
formalize the description of weak measurements. A path integral representation of 
continuous weak measurement described in this way leads directly to an analysis of 
backaction noise in terms of stochastic evolution equations. The average signal and the 
associated noise were obtained in terms of moments and time correlation functions of 
the measured quantity. 

In particular, the noise was shown to be an additive combination of a term 
characteristic of the measurement alone (detector noise) and terms associated with the 
system, which in turn include contributions from the intrinsic quantum mechanical 
uncertainty in the system and from backaction effects from the measurement process. A 
transparent representation of this stochastic evolution was obtained by separating it into 
a process characteristic only of the weak measurement, and another, representing the 
quantum uncertainty intrinsic to the system as well as that arising from the measurement 
backaction. This yields the noise as an additive combination of the corresponding 
contributions, while the total probability is found to be convolution of white Gaussian 
detections noise and intrinsic system's quasiprobability. The quasiprobability can be 
negative although the negativity is not visible at the level of second order correlations 
due to weak positivity. The general formalism was applied to two simple problems: 
continuous monitoring of the level population in a 2-level system and continuous 
measurement of the position of a harmonic oscillator. For both systems we have 
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established limits on the possibility to observe oscillatory motion in a single run of 
an experiment. The negativity property of the quasiprobability can be demonstrated on 
in the 2-level system, using fourth order correlations. In this way, we have constructed 
a Leggett-Garg-type inequality without the assumption of dichotomy or boundedness 
of the variable. 

We observe that the QZE occurs when both the Hamiltonian and the observable 
can be represented in finite-dimensional Hilbert space. When the space is infinite or 
continuous and both the Hamiltonian and the observable have no finite-dimensional 
representation, the dynamics will not always be able to 'pin down' the state and 
consequently the dynamics may get diffusive. Establishing criteria for the occurrence 
or absence of the QZE in realistic systems continues to be an intriguing and challenging 
issue. 
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Appendix A 

Here we derive ([5]). Using ([2]), a succession of time evolutions and measurements in the 
interval (0, t) reads 



K ({ 0i }) = (2A/ 7 r) Af /V (i/?l) ^ (tjv + 1 -* iv) e-^ a ^- i)2 ■ ■ ■ x 

e -(i/h)H(t 3 -t2) e -Ma2-A) 2 e -(i/h)H(t2~t 1 ) e -Ma 1 -A) 2 e -(i/h)Ht 1 



(A.l) 



Putting A = XAt and tj — tj-x = t\ = At and using At — > leads to 



k({ aj }) = (2AAt/ 7 r) JV / 4 f[e[- (i/r ^ A ^- i ) 2 ] A V^ 



(A.2) 




K[a{t)] = CTe^\r {l/r ^^ t) - A ) 2 \ 



(A.3) 



Alternatively, using A(tj) = eW** 1 * Ae~W dt ' yields 



k({ aj }) = (2XAt/n) N/i e- (i/h) 



N+l 




A(a J - J 4(t j )) 2 l At 



(A.4) 



and in the continuum limit 



K [a (t)} = Ce-WQ Ht Te- x fi( a ®- A ®) 2dt 



(A.5) 
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Appendix B 



Here we prove (10). Start from Texp J i(%(t)/2 + (fi(t))A(t)dt and discretize to get 



f e iJ(x{t)/2+<(>{i))A(t)dt _ j- e iAtE J (xfe)/2+0fe))A(i J ) = |^ = ■ _ ^ jyj 

e iAt( X (tN)/2+<i>(tN))A(t N ) e i&t(x(tN-l)/2+HtN-l))MtN-ri . . . e iAt(x(ii)/2+0(ti))A(ti) 
_ e (i/h)Ht Ne iAt( X (t N )/2+<l>(t N ))A e -(rmi/h)Ht N x 



Next replace 



e (i/ft)Ht N -i e iAt(x(*JV-i)/2+^(tAr-i))A e -(i/ft)m JV -i 
. . . e (i/R)Hti e iAt( X (ti)/2+^(ti))i e -(i/R)i?ti 



J iAt(xfe)/2+^.(t j )) J 4 _^ e iAtA x fe)/2 e iA 9 i(t J )(t,-^_ 1 ) 



for j = 1, ...N and define Ify (t) = H - hcj> (t) A, 

A^t) = f e m)So*6^f)dff •^ T - e -(i/R)/*^(f)df and discretize it again, 

A^fc) = f JJ e (i/ft)^(*,)A*i r JJ e -(i/fi)H ( tj )At 

to get 



(B.l) 
(B.2) 

(B.3) 
(B.4) 



from which (10) follows 



Appendix C 



Here we derive Eq.(19a 19c). Start from (13) and take its functional derivatives 



<a(t)> 



Da a (£) Pb[a] 



1 /oe Sf >M 



* V $X (t) / x(t)=0 
y D0e^ d ^ 2 «/ 2A Tr T (AS)f> + pM*)) / 2 

J D0e^ d ^ 2 «/ 2A Tr(i (t)p' 



(CI) 



which is (19a) 



(a (t) a (t')) b = Daa(t)a (t') P b [a] 



S 2 e s bix] 



$x(t)5 X (t')J x{t)=0 



\ J D^e^^Tr {A,(t'),{A,(t),p}} 



(C.2) 



Time ordering implies that for t' > t, t' will be placed in the outer commutator, 
however last expression is equal to 



i [ D0e-/^W/ 2A Tr \A^A^))p 



(C.3) 

which does not depend on the operator ordering. Higher moments are obtained in the 
same way. 
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Appendix D 

Here we consider the case [A, H] = and [p, A] = 0. When the observable A commutes 



with the Hamiltonian, A^ (t) = A so the trace in ( 13 ) becomes independent of 4>. Using 
/D0e-^ 2( * )/2A = 1 it follows that e s " w and P b [a , (pk, do not depend on A. This 



implies that the evolution associated with the backaction effect is deterministic and the 
only source of noise is the detector. To see the implication of this on the moments 



consider the moment generating function (c.f. (12), (13)) 

-/dt x 2 (t)/8A 



e s d [xl e s b [x] 



D0e" 



Jdt0 2 (t)/2A 



X 



Tr 7" e i /x(<)A (t)dt/2 p j- e i/ x (t)i (t)dt/2 
For the imposed initial conditions this becomes 



/dt x 2 W/8A (\af x (t)dt/2\ 



The first moment satisfies 
(a) 



,s d [x]+s b [ x ] 



x(t)=o 



% (t) 



AM 



a 



(DA) 



(D.2) 



(D.3) 



x(t)=o 



which implies also (a) = a when used in (22). It is easy to realize that the second 
moment satisfies 

5 2 



(a (t)a{t r )) : 
5 2 



s d [x] 



,S d [ X ]+S b [ X ] 



8 2 



x(t)=o 



M(t)6 X (t>y j x(t)=0 \5 X (t)5 X (t>y j x(t)=0 (D,4) 

The second term yields a 2 , so (6a (t) 5a (t')) = (5a (t) 5a (t')) — a 2 is determined just by 



AM 



the Gaussian detector noise that in the case of (22) results in 
((5a) 2 ) = l/(4Xt ). 

Appendix E 



(D.5) 



Here we derive the master equation (28). Following the steps that lead to ^ but 
without the trace, we find 



pw = e <^ m I D0e"^ 2 M/ 2A / D X e" / d ** 2 «/ 8 V /WW*)* x 



and 



j- e fi(x(t)/2+<f>(t))A(t)dt p f e fi( x (t)/2-<j>(t))A(t)dt e (i/h) 



p(t) = / Dap[a 



Ht 



(E.l) 



-{i/K)Ht 



D0 e " - 1 '^M/^j-jf <t>(t)A{t)dt p f e -if 



4>(t)A(t)dt Ji/h)Ht 



(E.2) 
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In what follows we will use the incremental propagation version of this equation: 

p(t + At) = e -mm+M) J D0e -/;+**d^(*')/2A x 

7V ft t+At Ht')Mt')dt ~ U\ f e -i P t +At Ht')Mt')dt e (i/h)H(t+At) 
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(E.3) 



Next use 



Te ijr^mMt')dt = Texp hAt<p(t)e^ m Ae-^ B 



(E.4) 



JJ e e/ft)£t Je iAt#t)A e -(i/R)ati = e (i/ro^(t+At) e / t t+At (i/r l )[0(t')i-H]dt' e -(i/s)^ 



to rewrite ( |E.3[ ) in the form 

^(t + At) = |D0e-^ +Atd ^ 2 W/ 2A x 
e i/ t *+ A «(^f)A-^)dt^^ e -i// +At (0(i')A-J/).!/ 



(E.5) 



We next expand the RHS of (E.5), keeping only terms that can contribute to order 
O (At) . To this end, we use 

e ±ifr*Ut(<p(t)A-H/H) = 1±L 4 [ t+At ^dtTifH/h) At 



t+At r t+At 



-A 2 / dtdt'(j){t)(j)(t')/2 

Jt Jt 

This leads, using (0) = and ((f)(t)4>(t' )) = A5(t - t'), to 

Pit + At) = pit) ~ *t[H, p{t)] (i/h) - i [A, Pit) 

-t+At 



(E.6) 



t+At 



(<Kt))dt 



+ 



(0(t)0(t')}(^(t)A-{A 2 ,p(*)}/2) 



pit) - At[H, pit)] (i/h) - XAt{A , p(t)}/2 + XAtApit)A 



(E.7) 



which yields 



dp(t) 



[H,p(t)\ - \{A 2 ,p(t)}/2 + \Ap(t)A 



dt h l 
= [H,P(t)]/ih-\[A, [A,p(t)}}/2. 



(E.8) 



Appendix F 



Here we derive (31 a 31c). We will demonstrate the derivation of the two-time correlation 



function, (|3 1 £>[) . We start from (j 1 9 £>[) in the form (|C.2|) and use the cyclic permutation 

(F.l) 



property of the trace together with identities such as 
to get 



(a(t)a(t')) b = - j I) A <p ! 



ATe-ft H <t>( s ) ds / ih x 



(F.2) 
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where H<f,{t) = H—h<f)(t)A. The functional integral J D\4> can be divided into a product 
of integrals performed over trajectories <j)(t) between and t and between t and t'. The 
former operates only on the dependent expression in the anticommutator brackets, 
yielding f D A 0fWo i ^( s ) ds / h pTe^o which satis fi es 



*(s)ds/H = (j ^ Q) p = ~ ^ 



(F.3) 



((F.3) is equivalent to (E.5), generalized for finite time evolution). (F.2) becomes 



(a(t)a(t'))i 



D A 0Tr 



Ate- 1 ^' M')*»/hjffi ( tj o) pTe^*' ^(»W fi l (F.4) 



Using (|F_3j again, now in the form j T> x ^Te~ ^ i6 ^ s)ds,h x {t)Te^ '^ {s)ds/h 
U (t',t)x(t), leads to (316). (31c) is verified analogously. 



Appendix G 



Here we prove (37). From (31a), (33) and (35), using also Ul = 1 (since LI = 0) and 
{frj, a k } = 25 jk l we get 

*(*), 



Tr 



^ (*))] 



-Tr (a 2 z ) =z(t). 



Next, from (316) and (33) 



((Tz(t)(Tz(t'))b 



4 



-Tr 



a 



,U(t',t){a z ,(l + a z (t))} 



(G.l) 



(G.2) 



for £' > t. Using (|35j) gives {(T 2 , (l + a z (t)) /2] =a z + z (t) 1. Finally, 

1. 



,U(t',t) (cr z + z(t) 1) 



{o- z {t)a z {t')) h = -Tr rr 
= ^Tr (a z (f -t) + z (t) 1)] i) 



(G.3) 
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